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Outer edges of face-on spiral galaxies * 
Deep optical imaging of NGC5923, UGC9837, and NGC5434 
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Abstract. We present deep optical imaging of three face-on disk galaxies together with a detailed description of 
the reduction and calibration methods used, in order to measure the intrinsic shape of their outer stellar edges. 
Whereas it is now well accepted that disks of spiral galaxies are not infinite exponential beyond galactocentric 
distances of about 3-5 radial scalelengths, the genuine structure of the truncation region is not yet well known. 
Our data quantitatively establish a smooth truncation behaviour of the radial surface brightness profiles and is 
best described by a two-slope model, characterised by an inner and outer exponential scalelength separated at a 
relatively well defined break radius. This result disagrees with the frequently assumed sharply truncated nature of 
the radial surface brightness profiles and implies the presence of stars and even star-formation beyond the break 
radius. In addition, we do not find a strong influence of a nearby companion on the ratio of the break radius to 
the radial scalelength. Our results denote new observational constraints for the search of the physical explanation 
for these smooth disk truncations. 
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1. Introduction 



Since the pioneering work of de Vaucouleurs (1948, 1958) 
and Freeman (1970) the radial surface brightness, tracing 
the stellar distribution, is empirically known to be well 
described by an exponential decline. The inner part can 
be dominated by a bulge, bar, or ring component, whereas 
the ou ter disk follows an exponential quite well with only 
minor deviations related to spiral arms. However, van der 
Kruit (1979|) has already pointed out that the profiles ap 
pear to break off towards the outer parts. This is later 
quantitatively confirmed using optical filters by e.g 



accepted there are still two major unsolved problems. 
On the one hand we are still missing a unique physical 
interpretation to describe this observational phenomenon. 
The proposed hypotheses span a rather wide range of pos- 
sibilities. Van der Kruit (1987) has deduced a connection 
to the galaxy formation process describing the trunca- 
tions as remnants from the early collapse. Ferguson & 



Clark (2001), for example, have proposed an evolutionary 



der Kruit fc Searle (19821 ), |Barteldrees fc Dettmar (1994| ) 



Pohlen ct al. (2000a|,H), and|de Grijs et al. (200l| ). Just re- 
cently, data by Florido et al. (2001) for the first time sug- 
gest these truncations to be also present in NIR images. 
These po called cut-off radii are practically only observed 
in edge -on galaxies due to the preferred liue-of-sighl inte- 
gration. While the existence of Lhe LmiicaLioiis is now well 



scenario represented by the viscous disk evolution mod- 
els. |Kcnnicutt (1989 ) has suggested a dynamical critical 
star- formation threshold, and Battaner et al. (2002) have 
proposed a connection with large scale galactic magnetic 
fields. 

On the other hand the detailed structure of the radial 
surface brightness distribution of disk galaxies is still not 



well constrained, although already Naslund & Jorsater 



(1997), Byun (199S), de Grijs et al. (2001), and Pohlen 



(2001) have found that the truncations are not that sharp 



Send offprint requests to: M. Pohlen, 
e-mail: pohlen@ll.iac.es 

* Based on observations obtained at the German-Spanish 
Astronomical Center (DSAZ), Calar Alto, jointly operated by 
the Max-Planck-Institut fur Astronomie Heidelberg and the 
Spanish National Commission for Astronomy. 



as thought before. While the forme r studies assume the 
truncation region to be exponential, Florido et al. (2001 ) 
used a model independent deprojection method to mea- 
sure and analyse the truncation curve (difference between 
measured profile and an infinitely exponential decline). 
However, nearly all studies used edge-on galaxies to search 
for truncations, whereas only the profiles of face-on galax- 
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ies are practically independent of possible line-of-sight 
effects caused by the integration across a multicompo- 
nent disk and less affected by the disturbing dust influ- 
ence. Observing face-on galaxies is therefore mandatory 
to study the detailed quantitative behaviour at the trun- 
cation region and to determine the scalelengths in an un- 
biased way. 

Up to now van der Kruit fc Shostak (1982| ) and shostak 



fc van der Kruit (1984) have reported the only known 



quantitative values of cut-offs in the literature for face- 



on gal axies (NGC 3938 and NGC628). |Van der Kruit fc 
Searle T982| ) have just found indications of a radial trun- 



cation for NGC 5033 and NGC 3198 and no evidence in 
NGC 4258 and NGC 5055, all being late type galaxies. 
Although calling the observed feature near the edge of 
the disk a 'truncation' van der Kruit fc Searle (1982j ) have 
already stated that "the edges are not infinitely sharp (we 
usually derive an upper limit to the scalelength beyond 
i? max of about lkpc)". In a subsequent paper van der 



Kruit 1 1988] ) h as stated that out o f the 20 face-on galax- 



ies observed by |Wevers ct al. (1986|) only four (NGC 2681, 
NGC 4203, NGC 5005, NGC 6340) —all being early type 
galaxies and therefore probably possessing a different disk 
formation history — do not show any sign for a drop off 
as judged from the last three contours. He has also recog- 
nized that only NGC 5371 shows a relatively sharp edge 
but without quoting detailed numbers. Explaining this in 
addition to the much lower surface brightness for face-on 
systems, with the problem of dealing with intrinsic devia- 
tions from the circular symmetry of disks which could be 
hidden by an azimuthally averaged profile. These asym- 
metries, such as spiral arms, especially from the young 
population should unavoidably contribute in less inclined 
systems. All these results are based on photographic mate- 
rial and up to now there is no systematic follow up of trun- 
cations for face-on galaxies. Especially a detailed quanti- 
tative modelling of radial surface brightness profiles based 
on higher quality CCD imaging data is missing. 

However, there are at least two recent papers in 
the literature which both have not found any sign for 
a truncation investigating two nearly face-on galaxies. 
Barton & Thompson (1997) — with a sensitivity limit of 



Mihn = 29.9mag/n" — have not found any evidence for 
a cut-off for the galaxy NGC 5383 out to approximately 
three Holmbcrg radii. They have measured an exponen- 
tially falling disk with several distinct departures from 
the smooth decline (two spiral arms with one outside the 
Holmberg radius), but their fitting yields unsatisfactory 
results. Nevertheless, in addition to a nearby compan- 
ion (UGC 8877) this galaxy unfortunately shows a rather 
complex morphological structure (type: PSBT3*P accord 
ing to the_RC30) and has also been rejected b y van der 
Kruit f l988 ) in his study. Wcincr ct al. (2001 ) also have 



found no sign of a truncation out to ten scalelength for the 




Fig. 1. Contour map (/x JR ) of NGC 5923 from 26.9 to 18.0 
equally spaced by 0.5 mag. 




1 Third Reference Ca talogue of bright galaxies (RC3: 
Vaucouleurs et al., 1991). 



Fig. 2. Contour map (/x JR ) of UGC 9837 from 26.8 to 18.0 
equally spaced by 0.5 mag. 



NGC 4123 (type: SBR5 + starburst according to NEDfJ) 
in their azimuthally averaged light profile. These studies 
have in common that both target galaxies are strongly 
barred and show a rather peculiar overall morphology. 

The detailed structure of the truncation region can 
be expected to give important constraints for the phys- 
ical nature of these outer stellar edges. In addition our 
analysis has consequential effects on the determination of 
the truncation as a structural parameter which is impor- 
tant to compare observational results with model predic- 
tions. Depending on the true physical nature these outer 
edges may contain fossil evidences imprinted by the de- 
tailed galaxy formation and evolutionary history. 

2. Sample Selection & Observations 

The main selection criteria for the observed objects are 
that they should be as close as possible to face-on disks. 
In a first step the galaxies are selected from the RC3 to be 
at least inclined less than 30° (deduced by cos(i) = a/b). 
According t o the lack of a truncation in the profile of S0- 
galaxies (cf. Pohlen, 2001 ) early type galaxies T < 2 are 



NASA/IPAC Extragalactic Database (NED). 
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(V 
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(9) 


NGC 5923 


152114.2 


+414333 


.SXS4.. 


4.0 


1.8 


5571 


5811.0 


80.7 


UGC 9837 


152351.7 


+580312 


.SXS5.. 


5.0 


1.8 


2657 


2934.0 


40.8 


NGC 5434 


140323.1 


+092653 


.SA.5.. 


5.0 


1.8 


4637 


4735.0 


65.8 



Table 1. Global parameters of the observed face-on galaxies: (l ) Principal name, (2) right ascension, and (s) declination, 



(4) RC3 coded Hubble-type, (5) the Hubble parameter T, and (6) the diameter in arcminutes ( de Vaucouleurs et al 
199l[)~[rhe (7) heliocentric radial velocities are taken from NED. According to the heliocentric radial velocities corrected 
for the Local Group infall onto Virgo (8) from LEDA, we estimated the (9) distances following the Hubble relation 



with the Hubble constant from the HST key project of Hq — 72 km s x Mpc 1 (Freedman et al., 2001) 
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Fig. 3. Contour map (/x JR ) of NGC 5434 from 27.0 to 18.0 
equally spaced by 0.5 mag. 

excluded as well as galaxies later than T > 6 to have a 
well defined and symmetric stellar disk. Additionally, we 
exclude galaxies classified as SB (including SABs) since a 
prominent bar may significantly disturb the assumed az- 
imuthal symmetry while smoothing the profiles. In order 
to obtain the selected galaxies with a decent resolution we 
first restrict the size to D25 2.0'. After visually inspect- 
ing all galaxies from the RC3 matching these criteria on 
DSS^| images, our target galaxies should possess as few as 
possible foreground stars within the FOV (especially no 
bright ones) and should not show significant irregularities 
pointing to present or recent interactions. It turns out that 
for a given observing period galaxies corresponding to all 
the above mentioned selection criteria and good visibility 
for approximately 4 h during the night are rare. Whereas 
from the visual inspection NGC 4136 seems to be the ideal 
galaxy, technical flatfieldin g pr oblems recognised at the 
telescope (described in Sect. |3.l|) forced us to choose galax- 
ies with smaller diameters and we relaxed the size limita- 
tion to D 25 > 1.8'. Therefore we finally choose NGC 5434, 
NGC 5923, and UGC 9837 as good candidates (cf. Fig.@, 
Fig. |, and Fig.|). NGC 5923 and UGC 9837 match all 
our criteria, despite some intermediate bright stars in the 
FOV, whereas NGC 5434 possesses two rather bright stars 



Digitized Sky Survey (DSS) 



within the FOV. Additionally, for NGC 5434 there is an al- 
most equal mass companion (NGC 5434B; KPG 410) with 
similar velocity nearby. Although NGC 5434 is therefore 
probably interacting and even slightly overlapping with 
the companion (in projection), its shape is extraordinary 
regular and free of large scale asymmetries. Owing to the 
lack of better suited galaxies, NGC 5434 was observed to 
investigate effects of a possible ongoing interaction on the 
radial surface brightness of both galaxies. The global pa- 
rameters of the observed face-on galaxies are shown in 
Tab.g. The images were obtained during an observing run 
in May 2001 at the 2.2 m telescope of the Calar Alto obser- 
vatory (Spain) equipped with CAFOS, providing a field- 
of-view of 16' using a SITe CCD with 15 /im pixel size 
(0.53 "pix -1 ). Showing a nearly rectangular filter char- 
acteristic combined with a significantly higher peak ef- 
ficiency, the so called Roser R (RR, #648/168) filter is 
favoured to increase the signal-to-noise using the same 
integration time and was used instead of the standard 
Johnson R (JR, #641/158) filter. The different filter char- 
acteristics are shown in Fig.|J. In order to account for 
the crucial flatfielding problem of deep surface photom- 
etry, images for all possible flatfield strategies were taken: 
dome-, twilight-, and nightsky-natfields. The nightsky- 
flatfields were acquired without loosing dark time by ob- 
serving additional edge-on galaxies (covering less area of 
the chip) and significant dithering of the rather small ob- 
jects (D25 «2') across the large FOV (016'). During each 
night several Landolt standard fields (Landolt, 1992) were 
observed to accomplish the photometric calibration. 

3. Image Reduction & Photometric Calibration 

The 2048x2048 pixel SITe CCD chip (RON: 5.3 e", GAIN: 
2.7e~/ADU) was windowed to improve the quite large 
read-out time (3min) of CAFOS. Approximately 400 pix- 
els were excluded on each side vertical to the readout 
axis to preserve the overscan region. Additional 10 bias 
frames were taken at the beginning and end of each night. 
It turned out that there is a minor light leak within the 
CAFOS camera and bias frames taken with the dome slit 
open are on average « 0.2 — 2.9ADUs brighter, depending 
on the actual brightness within the dome. Therefore all 
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Fig. 4. Roser R (RR, #648/168) and Johnson R (JR, 
#641/158) filter curves. 

bias frames are inspected and according to our records 
rejected. The following reduction procedures are per- 
formed with standard IRAF^ tasks. The approximately 
80 columns of the overscan region are averaged and an 
overscan vector is fitted with an 11th order polynomial, 
using a ±3<r rejection. This is removed row-by-row from 
each image by cedproc. A bias image is made for each night 
out of the remaining useful 10-15 overscan-subtracted bias 
frames with averaging and zt3cr rejection to remove cosmic 
rays by zerocombine. The resulting bias images still show 
offsets of w 4.7 ADUs together with a large scale gradi- 
ent and a large scale column and row pattern structure. 
This master bias is than subtracted from all flat and ob- 
ject frames. Dark frames taken during the day show that 
the dark current is equal or less than approximately 0.3 
ADUs per 600 s exposure with no evident additional two- 
dimensional structure. It was not possible to reduce the 
light passing through the closed dome-slit to zero during 
the day. It is therefore not possible to decide if the re- 
maining counts are either from real dark-current or from 

4 Image Reduction and Analysis Facility: 
http:/ /iraf.noao.edu/ 



stray light passing into the camera, consequently no dark 
correction is applied. 

3.1. Flatfielding 

One of the limiting conditions for deep surface photome- 
try is the accuracy of the flatfield. It turned out that both 
domc-flatfields and twilight-flatfields are able to remove 
the pixel-to-pixel variation but on cost of a significant 
large scale gradient remaining on the flatfielded science 
frames. Depending on the position of the dome slit to- 
wards the telescope/camera system during the time the 
twilight-flatfields were taken (90°difference between dusk 
and dawn), the residual large scale pattern on the flat- 
fielded images changes in the same direction. Therefore we 
have to assume that this structure is caused by stray light 
of the illuminated dome during the twilight-flatficlding 
as well as the dome-flatfielding procedure. An additional 
pupil stop used in the 4th night to minimize this stray- 
light did not solve the problem. Since we took approx- 
imately 20 science images each night with significantly 
dithering across the FOV we have tried to use a nightsky- 
flat median combined from all long-exposure science im- 
ages during the night after an additive median scaling and 
a 3<7 lower and 2er upper sigma-clipping with imcombine. 
It turned out that even with imcombine's highly sophis- 
ticated scaling and rejecting algorithms it is not possi- 
ble to remove all the structure in the combined nightsky- 
flat. There is still some residual pattern in the center of 
the image remaining. Taking all the 55 science images 
(^cxp > 600s) from the first three nights within the RR fil- 
ter together yields a decent result for the object residuals 
and a similar STDDEV compared to the dome- and sky- 
flats, but at the expense of some residual underlying large 
scale structure. Another possibility is to isolate the im- 
posed large scale gradient by using dome- or twilight-flats 
with the imsurRt task and divide the normalised surface-fit 
(third order) back into the image again. This works nicely 
for the saddle surface introduced by the twilight-flatfields, 
but not for the ring-like structure of the dome-flatfield. 
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A major drawback of this method is that one has to en- 
sure that the object — compared to the structure of the 
underlying surface — is excluded from the fit. Although it 
is possible for the face-on galaxies with the imsurfit op- 
tion 'invcirc' to adjust the position of the galaxy center 
for each fit as well as a radius, the task fails for the rather 
large edge-on galaxies chosen. The overall quality of the 
combined image relies in this case on the quality of the 12 
to 17 individual fits to the sky-flatfielded single images. 
In order to minimise the number of crucial fitting pro- 
cesses we finally used a third method. By combining all 
55 dome-flatfielded science images with imcombine using 
a multiplicative exposure time scaling, an additive median 
scaling, and a 3a lower and 2a upper sigma-clipping, the 
fairly constant (over all nights) dome-flatfield structure is 
isolated. After a median smoothing (15 pixel) to increase 
the S/N ratio and normalising of the image, this struc- 
ture is divided back into the dome-flatfield images to yield 
the final single images without depending on the imsurfit 
quality. 

3.2. Combining 

We test three different ways to combine the individual (up 
to 18) object frames. One combining method is a slightly 
modified version of the mosaic task from the ARNICAQ 
package, where we introduce the full capabilities of the 
IRAF imcombine-task instead of the predefined combine 
task using the provided centering routine. The mosaic-task 
allows only full pixel shifts in x- and y-direction with- 
out rotation or magnification. In addition, we use the cus- 
tom written IRAF script gacombine[] based on the avail- 
able IRAF centering task imcntr together with IRAFs 
geomap/geofcran-tasks to perform the actual matching of 
all the frames to one reference frame. We try to use the 
full transformation with subpixel-shifts, rotation and dis- 
tortion of the pixels and a more restricted approach using 
only subpixel-shifts. The object of comparison reveals that 
all three methods yields rather similar results, although 
the FWHM of selected stars in the gacombined image us- 
ing the full transformation is slightly (0.1") better. To ac- 
count for possible changes in the overall geometry during 
the rather large (> 3h) observing period the gacombine 
task with full transformation is chosen to create the final 
combined image. The imcombine task parameters are ad- 
justed to average the individual pixels after multiplicative 
exposure time scaling, additive median scaling, and a 3a- 
clipping using the cedclip algorithm. The additive median 
scaling of the original images is essential since variation in 
the sky level (mainly caused by different moon positions) 
during the night are of the order of 0.4 mag. The final 
combined images are shown in Fig.^| with high contrast 
to highlight the residual background structure. 



5 External IRAF add-on package, provided by L. Hunt, 
Arcetri (Italy). 

6 Could be provided upon request. 
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Fig. 6. Histogram of all pixel values in the star-masked 
image of NGC 5923 centered on the background peak. 

3.3. Sky subtraction 

The combined images are visually inspected and fore- 
ground stars and galaxies are masked by hand with the 
IRAF task imedit using at least three different radii de- 
pending on the object size. Accurate sky subtraction is 
crucial for the slope of the derived luminosity profiles espe- 
cially for the low luminosity outer part. The actual prob- 
lem is how and where to determine the single value which 
is called "the sky value". It should be a mean value for 
all the remaining background pixels on the image still in- 
cluding the galaxy, fainter foreground stars, and remain- 
ing background sources. There are several different pos- 
sibilities to determine such a characteristic value within 
the available data reduction packages IRAF and MIDASQ 
They are all aiming to determine either the median (mid- 
dle value) or the mode (the most probable value) of the 
pixel distribution. However, it turns out that neither the 
'MIDPT'-value (an estimate of the median of the pixel 
distribution) nor the 'MODE '-value determined with the 
IRAF imstat task^J over the whole image are efficient to 
determine a reasonable value for the background. Both 
overestimate the true background significantly. Therefore 
we tested three different methods to approximate the most 
common value — assumed to be the representative sky 
value — within the final image. 

1. A mathematical mode (the most frequently occur- 
ring) value, determined by eye as the local ex- 
tremum within the histogram of all pixel values 
plotted with imhistogram (cf. Fig.||). 

2. The 'MIDPT'-value determined with IRAF's 
imexam on several 5x5 pixel rectangular regions 
centered on the cursor position visually pointing to 
'empty' sky regions. 

7 Munich Image Data Analysis System: 
http://www.eso.org/projects/esomidas/ 

8 Note: Lowering the standard binwidth of 0.1 improves the 
quality. 
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Fig. 7. Sky subtraction problem: Azimuthally averaged 
radial profile of NGC 5923 using the 'empty' region sky- 
value (dashed line), the histogram method (solid line), and 
the exact median (dotted line). 



3. An exact median determined with the MIDAS task 
STATISTICS/IMAGE, which sorts all pixels of the 
image by their values. 

The problem becomes obvious comparing the different re- 
sults in the case of NGC 5923: 2903.0, 2901.4, and 2905.2, 
counts respectively. In Fig.fj] the resulting radial profiles 
are shown using these values. The profiles start to dif- 
fer below a surface brightness of fx « 27.0 JR-mag/n" , 
whereas the actual lcr-background fluctuations — usually 
used to determine the depth of the image — reach /j, ss 
29.0 JR-mag/n" . We finally decide to use the histogram 
method, since this produced the most straight profile be- 
yond /i « 27.0 JR-mag/n" without an additional curva- 
ture, continuing the undoubtedly linear profile from the 
break radius to this point. This prevents introducing an- 
other slope just at our limiting surface brightness. The 
uncertainty caused by the visual determination is shown 
in Fig.||. In addition to the mean sky the two associ- 
ated justifiable extremes are plotted. This discrepancy 
slightly affects the determined slope of the outer profile. 
In this case the profiles start to deviate significantly be- 
low /Lt « 27.2 JR-mag/n" which marks the final accuracy 
of our photometry. 

3.4. Photometric calibration 



Several Landolt standard fields ( |Landolt, 1992 ) have been 
observed each night following at least one over several air- 
masses to calculate the extinction coefficients. The cali- 
bration coefficients are determined according to following 
equations: 

jb -JB = C ,J B + C!,j B • A + c 2 ,jb • (JB - JR) 
jr - JR = c , jr. + Ci,jr • A + c 2 ,jr • (JB - JR) 

For the fourth night no reasonable solution is found and 
only a zero point could be determined. The applied coef- 
ficients are shown in Tab. The fields were observed in 




Fig. 8. Sky subtraction problem: Azimuthally averaged 
radial profile using different values obtained for UGC 9837 
with the histogram method (cf. text). 




Fig. 9. Contour map (hjb.) of NGC 5923 after transfor- 
mation into polar coordinates between 26.4 and 17.9 in 
0.5 mag steps. 



Johnson R and B as well as the Reiser R filter to quan- 
tify the difference between RR and JR. We do not find 
a systematic correlation of the RR-magnitudes and JR- 
magnitudes with star colour or airmass. Therefore, we as- 
sume in the following the mean difference of ffljR-mRR = 
0.35 ±0.02 mag. 

3.5. Polar transformation 

Using the first order symmetry of the chosen ideal face- 
on galaxies an azimuthal averaging within predefined seg- 
ments increases the signal-to-noise ratio and allows to fol- 
low a smooth profile down to the limits of the photom- 
etry. In order to independently pick up the size of these 
segments the whole image is transformed into polar coor- 
dinates using an already existing C-progranQ The polar 
image of NGC 5923 is shown in Fig.^. By combining the 
individual rows — split with the IRAF task imslice and 
apply the powerful rejecting and combining algorithms of 
imcombine — we are able to mimic an azimuthally smooth- 
ing within segments of the galaxy disk and to reduce the 
influence from individual H n regions. NIR images would 



Courtesy L. Schmidtobreick. 
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Date 


Band 


Zero point 


Airmass 


Colour 






Co 


coef. ci 


coef. C2 


Nl: 26.05.01 


JR 


-23.993 ±0.098 


0.178 ±0.035 


±0.030 ± 0.095 


N2: 27.05.01 


JR 


-23.973 ±0.134 


0.103 ±0.085 


±0.009 ± 0.095 


N3: 28.05.01 


JR 


-23.987 ±0.012 


0.096 ± 0.060 


±0.004 ±0.012 


N4: 29.05.01 


JR 


-23.675 ±0.302 






N5: 30.05.01 


JR 


-23.986 ±0.015 


0.103 ±0.008 


±0.003 ±0.015 



Table 2. Photometric calibration coefficients used for each night. 



be more suited to study the underlying mass distrubution, 
although currently probably infeasible to these depths. 
The optical passband used is therefore contaminated by 
star- formation in and outside the spiral arms. These in- 
trinsic deviations within the galaxy's structure from az- 
imuthal symmetry are obviously present (cf. individual 
profiles in Fig.[n]). However, our azimuthally smoothing 
by combining the individual sliced rows does not signifi- 
cantly change the inner or blur the outer slope it rather 
averages the individual slop es in a reasonable ma nner, in 
contrast to the statement of van der Kruit (1988 ). 



4. Analysis & Results 

4.1. Derived parameters 

We limit the following analysis of the radial profiles to 
the 360° averaged profiles, except for NGC 5434 where 
only the 180° profile, avoiding the edge-on companion, 
are taken into account. The final combined profiles to- 
gether with all four 90° segments obtained from the polar 
transformed image (see Sect. |3.5| ) are shown in Fig. |lO|. 
In all three cases we find a common behaviour. The ra- 
dial surface brightness profile is best described by a two- 
slope model, characterised by an inner and outer expo- 
nential scalelength separated at a rather sharp so called 
break radius Rbr- The inner approximately 8 — 10 arcsec 
are dominated by a bulge component followed by a rather 
well described first exponential decline of 3.3 mag within 
p»41"and a second also quite well described second expo- 
nential decline of 4.4 mag (down to the observing limit) 
within «27". Taking into account the estimated distance 
this translates into an inner approximately 2.6 kpc wide re- 
gion dominated by the bulge component followed by a first 
exponential decline of radially rather different sizes for the 
three galaxies (NGC 5434: Ai?«7.7kpc, UGC9837: Ai?« 
12.1 kpc, NGC 5923: AR « 18.1 kpc) and a second slope 
within AR w 7.8 kpc limited by the outermost measured 
point. The inner and outer slope are determined with an 
ID exponential fit of the form I(R) — Iq sxp(— R/h) to 
the combined radial surface brightness profile using the 
IRAF task tlinear. The inner boundary is defined by the 
bulge light and the fit is performed down to /^«28.5JR- 
mag/n" while the break radius is determined by eye. The 
resulting fits are shown in Fig.[ll]. The resulting values 
are shown in Tab.^, for the break radius i?br (4) in units 
of arcsec, kpc (with the estimated distance D (2)), and in 



units of the inner scalelength h m (5), the surface bright- 
ness jttbr (s)) at the position of the break radius, as well 
as the ratio of inner to outer scalelength hi n /h out (7)). 
The total sizes defined as the last measured point of their 
disks are: 30.5 kpc (NGC 5923), 16.0 kpc (UGC9837), and 
22.3 kpc (NGC 5434). Although being probably an inter- 
acting system the derived values for NGC 5434 are in no 
way exceptional compared to the other two galaxies. 

4.2. Comparison with the literature 

The observed smooth nature of the radial truncation is 
consistent with earlier work on edge-on galaxies. Already 
early photographic profiles fitted by van der Kruit & 
|Searle (198l| ) (e.g. their Fig. 9 for NGC 5907) show the 
smooth nature of the truncation and reveal that their 
sharply-truncated model explicitly overpredicts the lumi- 
nosity profile at the truncation region. A similar behaviour 
is found in the recent work of Pohlen ct al. (2000b|) where 
they fit the sharply-truncated model of van der Kruit & 
[gcarlc (1981 ) to CCD observations. 



The derived two-slope model consisting of two expo- 
nentials separated at a b reak radius is als o discussed by 
Naslund fe Jorsater (1997|) and |Byun (19981). By analysing 
a deep optical image of NGC 4565 Naslund fc Jorsater 



(1997 ) find a similar two-slope shape for their major axis 
profile (cf. their Fig. 2). They state that exponentials also 
provide a good solution at both sides beyond the break 
radius and that there is no clear hint in their data of any 
further change in the disk scalelength at surface bright- 
nesses lower than jiy = 28mag/n" . The same experience 
is made by Byun (1998 ) in the case of IC5249. He also 
discovers that the major axis of IC 5249 demonstrates a 
peculiar two-component structure with an exponential in- 
ner disk, where the change of surface brightness is slow 
and also an exponential outer disk with a steep change 
of surface brightness. He concludes that this structure 
cannot be caused by internal absorption but must rep- 
resent the intrinsic stellar distribution. Prior to this al- 
ready |Jensen fc Thuan (1982Q for NGC 4565 and [gasak 



(1987) for NGC 5907 fitted additional exponential profiles 



to the truncation regions. A different approach is outlined 



by Florido et al. (2001). They use a model independent 
deprojection method to measure and analyse the trunca- 
tion curve which is the difference between the measured 
profile and an infinitely exponential decline. 
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Table 3. Results from the ID exponential fitting of the radial surface brightness profiles using fixed inner and outer 
boundaries and tenably varying the break radius around the mean value. 




R [arcsec] R [aresec] 



Fig. 10. Combined 360°-segment (thick solid line) to- Fig. 11. One dimensional exponential fits to the two dis- 

gether with the four 90°-segments (thin dashed lines) tinct parts separately overplotted on the radial profiles: 

Radial profiles of NGC 5923 (upper panel), UGC 9837 NGC 5923 (upper panel), UGC 9837 (middle panel), and 

(middle panel), and NGC 5434 (only 180°) (lower panel). NGC 5434 (lower panel). 
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In the case of face-on galaxies van der Kruit & Shostak 
(1982|)]measure for NGC3938 a quite similar radial pro- 
file (cf. their Fig. 10) with a break radius at i?b r = 4.2-/i 



(fJ-br ~ [24.5 F-mag/ ) and a value for the ratio of inner 
and outer scalelength of h- m /h out — 1.6 which is consistent 
for the occurrence of the break radius and only somewhat 
smaller for the scalelength ratio compared to our results. 
Using the radial velocity of 1011 km s _1 from LEDAj^ cor- 
rected for the Local Group infall onto Virgo and a Hubble 
constant of Hq = 72kms~ 1 Mpc -1 (Frecdman ct al., 
200l|)7|we estimate the distance of NGC 3938 to 14.0 Mpc. 
Thereby we recompute the linear parameters and recover 
a rather low inner scalelength of 2.5 kpc, a more typical 
outer scalelength of 1.5 kpc, and a break radius at 10.2 kpc 
which are all similar to our values. In contrast to van der 



get galaxies. It is well known that the spiral structure is 
constrained between the inner Lindblad resonance (ILR) 
and the outer Lindblad resonance (OLR) (eg. CanzianJ 



1998). Unfortunately all three galaxies are not well de- 



Kruit 1 1 Shostak (1982 ) we do not regard our derived outer 
scalelengths as upper limits of the true scalelength. They 
assume this upper limits, due to the azimuthally aver- 
aged deviations from the pure circularity of the isophotes. 
It is obvious by comparing the 360° averaged profile to 
the individual 90° wide strips (cf. Fig. [l0|) that the az- 
imuthally averaging does not significantly change the mea- 
sured outer scalelength. There are of course variations be- 
tween the individual profiles due to spiral arms, prominent 
H II regions, and intrinsic asymmetries, but judging from 
the 360° averaged profile it seems to be a reasonable mean 
representation. 

The radial surface brightness profile of NGC 628 pre- 
sented by Shostak fc van der Kruit (1984[ ) also fits quite 
well with our results (cf. their Fig. 13). They find a break 
radius at precisely i?br = 3.1/ii n (/Xb r ~23.9F-mag/n" ), al- 
though quotuig_to_be_at^ 



fined grand-design spirals. NGC 5923 and NGC 5434 show 
nearly similar behaviour starting with a two-arm spiral 
which is splitting up towards the outer parts in at least 
four traceable arms. UGC 9837 is even more complex. In 
the inner part starting with a two-arm spiral, it looks 
rather flocculent in the outer parts. The estimated extent 
of the spiral arms are 60.4"(NGC 5923), 58.8" (UGC 9837), 
and 46.1"(NGC5434). In the case of NGC 5434 this fits 
well with the derived break radius, although for NGC 5923 
it is «6"(2.5kpc) an d for UGC 9837 » 10" (1.9 kpc) larger 
than i?b r • However, Canzian (1998 ) already points out 
that it is unknown how often spiral structure remains rel- 
atively bright, or even optically detectable, at radii ap- 
proaching the OLR. He derives for his sample of face-on 
to moderately inclined galaxies a large scatter of the spiral 
arm surface brightness at their outer extent, ranging from 
/i = 22.6 — 26.5 V-mag/n" , compared to our mean value 
of ^br =24.6R-mag/n" . 

Kennicutts study on the dependence of the star for- 
mation rate (SFR) on the density and dynamics of the 
interstellar gas (Kennicutt, 1989) proposes a dynamically 
critical gas density acting as a star formation thresh- 
old. Beyond this critical radius the star formation should 
be inhibited and if this persists for sufficient time, it 
should introduce a visible turnover or truncation of the 
observed stellar luminosity profile at that radius. These 
critical radii for star formation are measured by prominent 



the rat io of inner and outer scalelength of hj n /h out = 2.6, 
together with an inner scalelength of 4.0 kpc and an outer 
of » 1.5 kpc. 

Comparing our results with other published radial 
profiles of NGC 628 ( [Natali et al., 1992| ) also supports 
the two-slope model. According to our mean value of 
Rbx /hi B = 3.9 ± 0.7 and using their derived scalelength 
/iR-band = 3.3 we predict a break in the profile to be vis- 
ible at 13 ± 2 kpc. Assuming their profiles to be reliable 
out to 16 kpc (cf. their Fig. 8 and Fig. 9) we find a clear 
break in the B, V, R and I profiles with similar shape 
compa red to our three galaxies at 14 kpc equivalent to 



breaks in the Ha surface-brightness profiles (cf. Martin & 
Kennicutt, 2001| ). |Lclievre fc Roy (2000| ) discover that the 



Ha surface brightness profiles in the case of NGC 628 di 
minish monotonically as a function of galactocentric dis- 
tance, dropping dramatically at 20 kpc. This fall-off in 
Sr a near the edge of the optical dis k (defined by the 25th 
B-band isophote) is also observed by Ferguson et al. (1998 ) 



for NGC 628 and also NGC 6946. |Lelievre fc Roy (2000| ) 
conclude that the standard Schmidt law 'fails' at this ra- 
dius equivalent to a gas surface density of w 4 M Q pc~ 2 , 
although they find additional Hn regions out to 27 kpc, 
and an H I disk extending even much further 



R hl /h i _ 4.2. Therefore their results fit well to our re- (1998[ ) do not compare their results with two profiles 



However, Lelievre fc Roy (2000 ) and Ferguson et al 



suits for NGC 5923 and NGC 5434. Unfortunately their 
colour profiles just stop at this radial distance and we are 



tracing the underlying old stellar population which are 
present in the literature in the case of NGC 628. Shostak 



not ab le to find a similar break in the colour profiles which 



would possibly constrain the physical hypothesis for the 
origin of these truncations. 

4.3. Constraints for the origin of truncations 

In order to search for a possible correlation of the break 
radius and the presence of spiral arm structure we visu- 
ally determined the extant of the spiral arms fnr tlio tax= 



fc van der Kruit (1984|) locate a break radius for NGC 628 
at 12.4 kpc (= 12.9 kpc for D — 10.4 Mpc), whereas we 



deduce 14 kpc from [Natali et al. (1992[ ). 

Compared to the 'Ha -truncation' at R « 20 kpc of 
Lelievre & Roy (200C) optical breaks at 13 kpc or 14 kpc 



clearly rule out a one-to-one correlation. Nevertheless, at 
14 kpc a small plateau-like feature starts in the averaged 



Ha surface brightness (Fig. 3 of Lelievre & Roy, 2000). In 



addition, using the data for NGC 628 of Ferguson et al 



(1998) we deduce a clear break in the Ha profile at 0.8R 



2r, 



Lyon/Meudon Extragalactic Database (LED A). 



(cf. their Fig. 2). With their distance of 10.7 Mpc this is 
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equal to «13kpc (=12.6kpc for D = 10.4 Mpc) and there- 



fore similar to the break found by Shostak & van dcr Kruit 
(19840and the one deduced from Natali et al 



(1992) 



5. Conclusions 

Our deep optical imaging of three face-on galaxies con- 



firms the finding of Naslund & Jorsater (1997), Byun 
|(1998| ), |de Grijs et al. (2001] ), and |Pohlcn (2001| ) that the 
observed truncations of galactic stellar disks are not in- 
finitcsimally sharp but better described by an additional 
exponential decline. Since face-on scalelengths are not in- 
fluenced by line-of-sight integration problems as in the 
case for edge-on galaxies (so far investigated), we are able 
to quantitatively introduce a new two-slope model. We 
find that the radial surface brightness profiles are best 
represented by a two-slope exponential profile, charac- 
terised by an inner h{ n and outer /i out scalelength sepa- 
rated at a break radius i?t, r , rather than the previously 
assumed sharply-truncated exponential model proposed 



by van der Kruit & Searle (1981) and used for example in 
the study of |Pohlcn et al. (2000b] ). The mean value for the 
distance independent ratio of break radius to inner scale- 
length which marks the start of the truncation region is 
Rbr/h in =3.9 ±0.7. 

Although larger statistics are also definitely needed, 
the lack of any conspicuity in the parameters for 
NGC 5434 with its close physical companion points 
against an interaction hypothesis for explaining stellar 
truncations. In addition, we could not correlate undoubt- 
edly the end of the visual spiral arms with the occurrence 
of the break radii, which may be a hint relating the trun- 
cation phenomenon to star formation properties. Similar 
deep imaging of grand-design spirals are needed to clar- 
ify this hypothesis. Since we do not have deep Ha maps 
for our galaxies we have searched the literature for signs 
of break radii in published surface brightness profiles and 
have compared them in the case of NGC 628 with already 
available Ha data. Although not conclusive in this case an 
Ha follow up of the three observed face-on galaxies will 
be a promising project to determine the origin for the ob- 
served break radii (or smooth truncations) and address the 
star formation threshold hypothesis by Kennicutt (1989| ) . 

However, the detailed quantitative structure of the 
smoothness of the observed truncations presented here will 
set strong constraints for present and future theoretical 
models trying to explain this observational phenomenon. 
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